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ABSTRACT 

Context. Astrometric observations of resolved binaries provide estimates of orbital periods and will eventually lead to measurement 
of dynamical masses. Only a few very low mass star and brown dwarf masses have been measured to date, and the mass-luminosity 
relation still needs to be calibrated. 

Aims. We have monitored 14 very low mass multiple systems for several years to confirm their multiplicity and, for those with a short 
period, derive accurate orbital parameters and dynamical mass estimates. 

Methods. We have used high spatial resolution images obtained at the Paranal, Lick and HST observatories to obtain astrometric 
and photometric measurements of the multiple systems at several epochs. The targets have periods ranging from 5 to 200 years, and 
spectral types in the range M7.5 - T5.5. 

Results. All of our 14 multiple systems are confirmed as common proper motion pairs. One system (2MASSW J0920122+35I742) 
is not resolved in our new images, probably because the discovery images were taken near maximum elongation. Six systems have 
periods short enough to allow dynamical mass measurements within the next 15 to 20 years. We estimate that only 8% of the ultracool 
dwarfs in the solar neighborhood are binaries with separations large enough to be resolved, and yet periods short enough to derive 
astrometric orbital fits over a reasonable time frame with current instrumentation. A survey that doubles the number of ultracool 
dwarfs observed with high angular resolution is called for to discover enough binaries for a first attempt to derive the mass-luminosity 
relationship for very low-mass stars and brown dwarfs. 

Key words. - stars: low-mass, brown dwarfs, imaging - Binary: visual - Astrometry 



1. Introduction 

Over the last few years, intensive computational and observational efforts have been made to improve our understanding of the 
formation processes and evolution of brown dwarfs (BDs) and very low mass (VLM) stars. The determination of their Initial Mass 
Function (IMF) is a crucial step in this direction. Translating an observed luminosity function into an IMF requires an accurate de- 
termination of their mass-luminosity relationship at different ages, which up to now relies primarily on theoretical mass-luminosity 
relationships. Although the empirical constraints on these relat i onships for VLM stars have considerably improved within the past 
years (see e.g Hillenbrand & White, '2004; 'Delfosse et al.','200(y; Segransan et al.,'2000') only a few observational constraints a re cur- 
rent ly a vailab le and large uncertainties remain (Leinert et al., 2001; Bouy et al.. .2004b; .Brandner et al... 2004; Zapatero Osorio et al.L 
120041; IClose et al.l 120051: IStassun et al.Ll2006Hlreland et al.Ll2008l) . 
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The degeneracy in the mass-luminosity relation for ultra-cool dwarfs (UCDs) makes it difficult to accurately estimate their 
physical properties. Dynamical masses, which are not model-dependent, are a unique way to calibrate this relation. The components 
of a multiple system are expected to be coeval, removing part of the above mentioned degeneracy. Although the ages of the targets 
studied in this work are not well constrained, it will be possible, once their dynamical masses are known, to take advantage of 
their coevality to test the evolutionary models. By adjusting the theoretical isochrones empirically to fit both the observed total 
masses and the individual luminosities of the multiple systems, it will be possible to directly check the consistency of the models 
with the observations. The corresponding predictions on the age can then be compared to other indicators such as the activity, the 
rotation, and the presence and strength of particular spectral fea tures (such as L i, Ha), b ut also to more recent techniques based 
on spe c tral analysis of gravity sensit ive f eatures as described by iMohantv et all ( 2004). McGovern et alj (l2004t) . iMartm & Osorid 
(l2003h . IZapatero Osorio et all (|2004|) and iBurgasser et all (|2006|) . Finally, by studying the physical characteristics of objects with 
known dynamical masses, it will be possible to provide crucial information for our understanding of their physical properties, such 
as their interior structure, the formation of dust, the settling and depletion of refractory elements, and the underlying opacities. An 
accurate determination of the mass of an object based on dynamical masses in binary systems therefore provides not only a reality 
check for the theory but also a cornerstone in the understanding of the mass distribution of brown dwarfs. 

In this work, we present a time-series of high angular resolution observations aimed at monitoring binary ultra-cool dwarfs. 
These observations confirm the common proper motion of the binary candidates and represent a first step towards the derivation 
of orbital parameters and dynamical masses. Most objects presented here were monitored over timescales too short in comparison 
with their periods, allowing us to estimate rough orbital periods, but preventing us from obtaining detailed orbital fits. 



2. Observations and data analysis 

We have monitored 14 multiple systems using a variety of instruments; in the Hubble Space Telescope (HST), we used the Wide 
Field Planetary Camera (WFPC2), the Space Telescope Imaging Spectrograph (STIS), the Near-Infrared Camera and Multi-Object 
Spectrometer (NICMOS) and the Advanced Camera for Surveys (ACS). From the ground, we collected observations using the 
Very Large Telescope (VLT) with the NACO adaptive optics system, and the Lick Observatory Shane 3-m telescope, also with an 
adaptive optics system. Table[3]gives an overview of the characteristics of these instruments. The observations reported here started 
in 2000, and continued until the end of 2007. Table |4] lists the observations recorded per target, together with the corresponding 
relative astrometry, and photometry analysis. 

2.1. Sample 

The sample includes 14 binaries ranging from spectral M7.5 to T5.5 (see Table |2|i, therefore cove ring a wide range of primary 
masses. All objects wer e known to be multiple and had been observed and reso l ved at least once (see lReid et al.Ll2001tlBouv et al.L 
l2003HGizis etail 120031: IcTose et all 120031: IBurgasser etalL 120031: IsTederet al.Ll2005h . 

2.2. HST/WFPC2 observations 

We used the HST/WFPC2 and its Planetary Camera (PC, lBirettall2002l) with th e F814W filter (progra ms GO-9157,GO-9345 ,GO- 
9499 and GO-9968, PI. E. L. Martin). Part of these datasets were pubUshed in lBrandner et alj (|2004|) and lBouv etal] (l2005h . All 
targets were centered in the Planetary Camera (PC) which provides the best sampling of the PSF In order to identify and remove 
cosmic ray events and bad pixels, we used a four-point dithering pattern with typical exposure times of 400 s, adding up to a total of 
1600 s in total for each target. The images have been processed following standard procedures using the STScI STSDAS package 
in IRAF together with the calibration files provided by the STScI team. 

2.3. HST/NICMOS observations 

We retrieved NICMOS data from the HST public archive (program 9843, PI. Gizis). Two objects (2MASSW J0850359-I-105715 
and 2MASS W J 1728 1 14-H394859) have been observed using the NICMOS 1 camera. The objects were observed in MULTIACCUM 
mode with exposure times of 128 s and 144 s respectively. The data were processed following the recommendations of the HST 
Data Handbook with the STSDAS pipeline in IRAF and standard STScI calibration files. 

2.4. HST/ACS observations 

We started using the HST/ACS and its High Resolution Channel camera (HRC, iPavlovskvi l2003h with the F625W, F775W and 
F850LP filters (program GO-9451, P.I. Brandner). Later observations were collected using only the F814W filter (GO -10559. P.I. 
Bouy) t o obtain additional epochs for known VLM binaries. Part of the earlier datasets have already been published in lBouv et al.l 
(I2004allbh . The data were obtained in CR-SPLIT mode with a four points dithering pattern in each filter, and typical exposure times 
of 490 s, 230 s, 180 s and 300 s with the F625W, F775W, F850LP and F814W filters respectively. The images have been processed 
following standard procedures using the STScI STSDAS package in IRAF together with the calibration files provided by the STScI 
team. 
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2.5. HST/STIS observations 

As part of program GO-9451 (P.I. Brandner), spatially re solved STIS spectra of binaries of this sample were obtained using the 
high spatial resolution STIS spectrograph on-board HST jKim Ouiianoi l2063h . The corresponding spectroscopic data have been 
described in detail in another paper (Martin etal., 200^. This paper focuses only on the pre-acquisition images obtained with 
STIS prior to each spectroscopic exposure. These images were obtained in the Longpass filter (A^en =7230A, FWHM=2720A), 
with typical exposure times of 5 to 10 s. They have been processed following standard procedures as described in the STIS User's 
Handbook using the STScI STSDAS package in IRAF together with the calibration files provided by the STScI team. 

2.6. VLT/NACO observations 

We used the adaptive optics system NACO in order to obtain high spatial resolution images of VLM binaries (programs 70.D-0773, 
077.C-0062, 71.C-0327, P.I. Bouy). NACO and its near-infrared wavefront sensor provided excellent diffraction limited images of 
the binaries. Prior to period 71, we requested to use the AO system with the N20C80 dichroic. This dichroic allows 80% of the 
near-infrared light to reach the NIR wavefront sensor and 20% to be collected by the ALADDIN detector of the science camera. 
After period 71, all images have been obtained with the N90C10 dichroic that sends 90% of the light to the adaptive optics and 10% 
to the ALADDIN detector, allowing to close the loop on even fainter objects. Our scientific targets were used as reference star for 
the wavefront sensing. The images were obtai ned in jitter mod e with a four or five points dithering pattern. We processed the data 
with the recommended Eclipse jitter package ( lDeviIlardlll997h and the calibration files provided by the Paranal observatory. 

2.7. Lick/AO observations 

We used the adaptive optics facility of the Lick Observatory Shane 3 m telescope jGavel et al.L '2002") on 2007 April 4* to observe 
2MASS J 1847034-1-552243 (using H and Ks broad-band filters) and 2MASS J 1047 138-1-402649 (Ks only). These two targets and 
their neighboring stars are too faint to be used as natural guide star for Lick/AO wavefront sensing. Thus, we used the AO together 
with the Laser Guide Star (LGS) system. The Lick LGS system can perform tip-tilt wavefront sensing on a reference star brighter 
than R<16.5 mag. For the tip-tilt wavefront sensing, we used USNO-Bl.O 1304-021 1669 (ci'=10h47minl2.61s, (5=-)-40°26min44.0s, 
R=16.5 mag) and USNO-Bl.O 1453-027661 1 (a=18h47min00.6s, (5=+55h22min25.3s, R=15.6 mag), located at 14'.'4 and 29'.'8 of 
2MASS J 1047 138 -1-402649 and 2MASS J 1847034-1-552243, respectively. The laser spot was used for higher order corrections. The 
targets were observed using a 5 point dithering pattern, with exposure times of 30 s at each position. A PSF reference star was 
obtained just after 2MASS J1847034-H552243. In the case of 2MASS J1047138+402649, we used the first component of the system 
as reference PSF, ensuring optimized results for the PSF fitting procedure. 

2.8. Analysis of tine data 

In ord er to obtain the precise relative astrometry of these multiple systems, we used the same software described by iBouv et aTl 
(120031) . adapted to ACS, STIS, NICMOS, Lick/AO and VLT/NACO. The program, its performances and limitations are fully de- 
scribed in the paper cited above. A single point source can be described by only three parameters: the position of its centroid (x,y), 
and its total flux (/). A binary system is described by 6 parameters. The custom made program makes a non-linear fit of the binary 
system, fitting both components simultaneously rather than individu ally. It uses a library of 10 reference PSF (9 natural PSF and 
1 Tiny TIM synthetic PSF in the case of HST, iKrist & HooH |2003|) . except in the case of Lick/ AO and NACO, for which only 
one reference PSF star obtained the same night with the same instrumental settings was used. A ;t'^-minimization between the 
synthetic binary and the ob s erve d binary gives the b est values for the six parameters. Typical uncertainties and systematic errors 
are described in lBouv et"aLl (1200 3) and Bouvl ( l2004l) for both ACS and WFPC2. Similar calibrations have been done for Lick/AO 
and NACO. Briefly, for well resolved multiple systems with moderate differences of magnitude, uncertainties and systematic errors 
are in general estimated to add up to a; 10% of the plate-scale of the instrument, provided that 3 conditions are met: a) the PSF is 
well sampled, b) the reference PSF is of good quality and c) that the signal-to-noise ratio is l arge enough . The efi'ective resolution 
also depends on th e technique used to measure the relative astrometry and photometry (see e.g lClose et al.Ll2002tlBouv et al.L 120031: 
iKraus et all |2005, for 3 independent techniques). The values quoted in Table [3] are only indicative and relatively conservative. 
These uncertainties do not include systematic instrumental errors, which are discussed below and can sometimes dominate. Table [1] 
gives an overview of these systematic errors. They should be added quadratically to the uncertainties given in Table|4] 

2.8.1 . HST/WFPC2 systematic errors 

The main systematic errors on relative astrometry are due to: 

- the uncertainty on the absolute roll angle of the spacecraft (<0.003° according to the User's manual) 

- 34* row defect producing an astrometric offset of approximately 3% of the pixel height every 34 rows 

- the geometric distorsion (<0'.'005 of error according to the User's manual) 

The separations of the multiple systems presented in this paper are all less than 13 rows, so that the 34* row defect affects them 
once at most. The maximum systematic errors on the relative astrometry measured with WFPC2 therefore adds up to 0'.'0052, and 
the position angle to 0.003°. 
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Table 1 . Estimates of the maximum systematic astrometric errors obtained with HST, VLT/NACO and Lick/ AO 



lllall UillClll 


CllUl vjcp. 


n 


HST/WFPC2 


5.2 mas 


0.003 


HST/ACS 


2.8 mas 


0.003 


HST/STIS 


5.1 mas 


0.003 


HST/NICMOS 


4.4 mas 


0.003 


VLT/NACO 


1% 


0.31 


Lick/AO 


1% 


0.35 



2.8.2. HST/ACS systematic errors 

The systematic errors are primarily due to the accuracy of the roll angle of the spacecraft (< 0.003° as above) an d to the accuracy 
with which the geometric distortion of the camera has been characterized. The MultiDrizzle (iKoekemoerl l2005h pipeline corrects 
for most of the geometric distorsions, and the final relative astrometry is expected to be better than 0.1 pixel, or w0'.'0028. 

2.8.3. HST/STIS systematic errors 

As in the case of WFPC2 and ACS, the systematic errors are mainly due to the accuracy of the orientation of the spacecraft (<0.003° 
as above) and to the stability of optical distortion. The STI S Instrument Handb ook gives an accuracy for relative astrometry within 
an image better than 0.1 pixel, corresponding to xiS.l mas dKim Ouiianoll2003l) . 

2.8.4. HST/NICIV10S1 systematic errors 

The NICMOS pixel scales along the X and Y axes of each camera are slightly different, because of the slight tilt of the NICMOS 
arrays relative to the focal plane. The difference is of the order of 3%c only, and we neglect it in our analysis. The distortion 
corrections for the NICMOS 1 camera are small, even at the edge of the camera (0.9 pixels). After correction using the drizzle 
package provided by the STSci team, the relative astrometry in the center of the camera where all our targets were observed is 
expected to be better than 0.1 pixel corresponding to 4.4 mas. As in the case of the other HST instruments, the systematic errors 
also include the uncertainty on the orientation of the spacecraft (<0.003° as above .Barker. 2007.) . 

2.8.5. VLT/NACO and Lick/ AO systematic errors 

In addition to static instrumental uncertainties, images obtained with AO are known to suffer from variable effects, due in particular 
to temporal and spatial variability of the atmospheric conditions. These effects can vary significantly on short timescales, even 
between two consecutive exposures, and thus a recorded PSF is only an approximation of the system's PSF. In the case of Lick/ AO, 
we measured the effective platescale and position angle using a set of astrometric calibrators. The platescale was found to vary by 
as much as 1%, corresponding to 0.8 mas/pixel, and the position angle to be off by as much as 0.34°. Our NACO observations were 
made in service mode with stan dard calibrations, and no astrometric calibrators were therefore obtained to control the platescale 
and orientation accuracv. Eggenberger et al] ( l2007h report recent measurements of the instrumental uncertainties obtained for NACO 
with similar settings. They measure platescale variations as large as 1%, corresponding to 0.14 mas/pixel, and position angle offsets 
as large as 0.31°. 



3. Analysis 

3.1. Common proper motion pairs 

Ten objects have proper m otion measurements in I Jameson et al. I (l2007h . lDahn et al.1 (l2002h . lTinnev et alJ (l2003h . lVrba etaP (l2004 
or the USNO-B.l catalog (iMonet et al.L l2003h . All but three of these targets are confirmed as common proper motion pairs with 
motion of the secondary much lower than the proper motion (see Table |5]l. For the L-dwarf pairs 2MASSW 2331016-040619 
and 2MASSW J1728114-I-394859, Table |5] gives a proper motion amplitude comparable to the motion of the secondary, but the 
orientation of the proper motion of the unresolved pairs is inconsistent with the companion being an unrelated background source, 
as illustrated in Fig.H The L dwarf 2MASSW J0920122-H351742 is not resolved in our new VLT and HST images (see Section[l2l 
for a detailed discussion on that particular object). Even though accurate kinematics measurements are required to confirm that the 
objects without proper motion measurements are comoving, we note that the motion of the secondary component with respect to the 
primary is consistent with that expected for a gravitationally-bound companion. Considering the uncertainties, the relative motion 
is of th e order of a^lO mas/yr, typically lower than the proper motions expected for such nearby objects (^100 mas/yr, see Table|5] 
and e.g lDahn et aLll2002tlTinnev et al.Ll2003h . 
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Fig. 1 Mosaic of images of 2MASSW J0920122-I-351742. The observation date and instrument are indicated. The scale is repre- 
sented and is the same in each image stamp. Contour plots are over-plotted to illustrate the clear elongation in the first epoch image 
(.Reid et al..,2001.) . the possible elongation in the last epoch image in the same direction, and the round PSF at the other epochs. 



3.2. 2MASSWJ09201 22+351 742 

2MASSW J0920122+351742 (L6.5) has been unambiguously resolved as a binary by lReid et all (1200 ll) using HSTAVFPC2, with 
a separation of 0'.'075. This object is not resolved by us, neither in our 2 epochs with HST/ACS and HST/STIS, nor in our third 
VLT/NACO epoch. Figure[T2lshows a mosaic of the 6 epoch images of 2MASSW J0920122-f35 1742 obtained with HST and VLT. 

The object is clearly elongated in the WFPC2 image, as shown in Figure [1] It is elongated in the three consecutive images 
obtained that day, excluding the possibility of a cosmic ray event. Moreover, other objects present in the field of view of the WFPC2 
images do not show any elongation, excluding any instrumental problem. 

The presence of the nearby stai- 2MASS J09201092-f35 17452 in the Februai-y 2000 WFPC2 image and in the March 1998 
2MASS images allows us to rule out the combination of a high proper motion brown dwarf with a background star aligned by 
chance at the first epoch. 

Figure [U shows that the PSF of the 4 consecutive ACS (2002), STIS (2003), NACO (2003) and ACS (2005) images look sharp 
and unresolved. Because the system is not resolved, we can put an upper limit of ^0'.'06 on the separation of the two components 
of the system, corresponding over the 5.6 yr time difference to a motion of O'.'Oll yr ~^ This measured motion is much sm aller 
than the typical O'.'lOO yr ' reported for such nearby ultracool dwarfs (see Table |5] and iDahn et all 120021; iTinnev et al.L 12003) . and 
suggests that the absence of motion detection is due to the fact that the pair is comoving (assuming negUgible motion for eventual 
background coincidence). An accurate proper motion measurement should confirm this preliminary conclusion. 

A more detailed analysis of the last epoch image (2006) obtained with ACS shows that the PSF seems a little elongated. 
Figure |2] shows a comparison of the residuals after PSF subtraction of the resolved WFPC2 image, the unresolved ACS image of 
2005 and the possibly resolved last epoch ACS image. The residuals are significantly stronger in the first and last one, with an 
elongation in the same direction, indicating that the object is possibly al most resolved in the last epoch. The first and last epochs 
are separated by 5.6 yr, close to the estimated orbital period (a:7.2 years. IBouv et all '2003). This suggests as possible explanation 
that the companion might have been too close to be resolved in the NACO, ACS (2003 and 2005) and STIS images, while close 
to its maximum elongation in the WFPC2 and last ACS image. The relatively short estimated period of a!7.2 years, and the short 
separation (only 0'.'075, very close to the limit of resolution of HST and VLT at these wavelengths) are consistent with such a 
scenario. 

Simple calculations considering an eccentric orbit, with a period of 7.2 yrs, a semi-major axis of 0'.'075, as measured in the 
WPFC2 image, and the companion at its apastron at the date the WFPC2 images indicate that the probability that the companion 
could not be resolved by either NACO, STIS or ACS is relatively high. Figure[3]illustrates these calculations in the cases of typical 
eccentricities of 0.1, 0.3 and 0.5. In these configurations, and for eccentricities greater than 0.3, the companion would have been 
resolved (or almost resolved) in the last ACS image but in none of the other ACS, STIS or NACO images. Although simplistic, 
these calculations show that further observations of 2MASSW J0920122-I-351742 will have to be taken near maximum elongation 
in order to resolve the binary again with currently available instruments. 
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Fig. 2 Comparison of the average residuals obtained after single-star PSF subtraction on the resolved WFPC2 image of 
2MASSW J0920122-I-351742 (left), the unresolved ACS image (middle), and the possibly resolved ACS image (right). The color 
scale and orientation (North/Up and East/Left) are the same in each image. The scale is indicated in the left stamp and is the same 
for each image. The residuals are significantly stronger in the 2000 WFPC2 and in the 2006 ACS image than in the 2005 ACS 
image. 
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Fig. 3 Figure illustrating possible scenarios explaining why the companion was detected in the HST/WFPC2 images but 
not in the VLT/NACO, HST/STIS and HST/ACS ones. It assumes orbits with typical eccentricities of 0.1, 0.3 and 0.5 for 
2MASSW J0920122-H351742AB, with a semi-major axis of a.'075, a period of 7.2 yrs, and the apastron at the epoch of the WFPC2 
observations. The filled circle in the center indicates the position of the primary, and the diamonds the position of the secondary 
at the epochs of the 5 observations. The horizontal dashed and dash ed-dotted lin e s indicate th e limit of resolution of respectively 
WFPC2/STIS (~0'.'060), and ACS/NACO (~a.'045), as estimated in lBouvl (l2004h : iBouvetan (i2003.) . The line of sight, chosen in 
the most favorable case, is indicated on the right. 



3.3. Discussion on particular objects 

2MASSW J233 1016-040619 - Table |4] and Fig. |4] show that the consecutive measurements obtained for this multiple system do 
not follow a keplerian motion. The Gemini Hokupa'a measurement taken in 2001 indicates a clockwise motion, while the VLT 
NACO 2003 measurement suggest a counter-clockwise motion. The NACO images, with an exposure time of only 2x30 s during 
very poor ambient conditions, were of low quality compared to the Gemini and HST images. The faint companion (AK=2.44 mag), 
is barely detected in the NACO images, and the corresponding astrometric measurement is therefore not reliable. We report it for 
completeness, but it should be considered with caution. 

DENIS-P J035726.9-441730- Tableland Fig.|4]show that the separation was increasing until 2003, and decreasing in the following 
years. If confirmed (the uncertainties are relatively large), this would mean that the observations were obtained close to the maximum 
elongation. Assuming a nearly edge-on orbit, as suggested by the current measurements (see Fig.|4]i, a photometric distance corrected 
for multiplicity of 22.2 pc (Bouv et al., 2003), and a total mass of 0.15 M^^j, the semi-major axis of » 104 mas corresponds to a period 
of ~9 yr If the distance is larger, i.e., 30 pc, the period remains relatively short (a!l5 yr). Hence, DENIS-P J035726. 9-441730 is a 
promising target for dynamical mass measurement within the next few years. 
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3.3.1. Period estimates 

Table |6] gives estimates of the orbital periods using three different methods. The first column gives the p eriod calculated using 
Kepler's third law, a statistical scaling of the measured angular separation by 1 .26 (! Fischer & Marcvi 1199 2'). trigonometric, photo- 
m etric or spectra-ph otometric derived distances, and mass estimates derive d from the spectral types using the spectral type vs Teff 
of iDahn et al.l ( |2002|) and the evolutionnary models of iBaraffe et al.l (1 1998b for an age of 1 Gyr. The second column gives periods 
calculated from the fractional change in P.A assuming a circular face-on orbit. The third column gives the period calculated using 
the ratio of the total motion over the monitored timespan to the quantity 4xmaximum projected separation, assuming an edge-on 
circular orbit (an edge-on circular orbit would traverse the maximum separation approximately four times). Although the hypothesis 
are strong and numerous and the uncertainties large, the agreement between the different estimates can sometimes provide an idea 
of the real value of the period, as well as some idea of the inclination or eccentricity of the system. A large difference between 
the second and third column can indeed imply a large inclination and/or eccentricity. The cases of 2MASSW J15344984-2952274 
and 2MASSW J1426316+1557013 illustrate the effect of inclination. For these two pairs, the second column (face-on case) gives 
a period estimate respectively ^12 and 50 times larger than the third (edge-on case) and first columns. A comparison with Fig. |4] 
shows that these systems are seen nearly edge-on, explaning the strong discrepancy. In some cases, the discrepency between the 
first column and the two other gives some clue on the eccentricity. In the case of DENIS-P J035726. 9-441730, for which Fig. |4] 
shows that it was observed at the apastron passage, the large discrepancy between the first column and the other two indicates that 
the system most probably has an eccentric orbit. The effect of inclination and eccentricity being degenerated, and the uncertainties 
on the distances and masses being large, it is not possible to go beyond such qualitative discussions. 

4. Future Prospects and conclusions 

We present astrometric and photometric results of follow-up observations of 14 UCD binaries. Only half of them are rotating fast 
enough to provide accurate dynamical masses within the next 15-20 years. The HST, but also the recently commissioned Laser Guide 
Stars for Adaptive Optics on 8 m class telescopes should allow to discover and follow more UCD binaries, usually too faint and too 
red even for the IR-WFS of NACO. Some targets not included in the present sample are already part of other on-going programs, 
and more follow-up observations are likely to be published in the coming months/years. We are currently closely monitoring three 
additional targets for whic h dynamical masses will be derived within one year ( Bouy et al., in prep.). Another two (e-Indi Bab and 
GJ lOOlBC, respectively iMcCaughrean et aTl 120041; iGolimowski et al.l |2007[) are the targets of additional monitoring programs. 
The total number of "short" period VLM multiple systems (short meaning periods allowing dynamical mass measurements within 
15-20 yr) ro ughly adds up to a dozen of obj ects, which has been extracted from original samples of UCDs made of xl40 objects 
jBouv et all 120031: IClose et al.l 120031 : iGizis et al.l l2003b . i.e. the frequency of short-period resolved binaries is about 8%. If we 
consider that about 20 binaries (40 masses) are required in order to start calibrating the mass-luminosity relationship, the current 
study shows that we would need to observe a total of roughly 140/12x20=250 UCDs at high spatial resolution. This estimate means 
that another survey of about 140 more UCDs is needed to discover enough binaries that can yield dynamical masses in the near 
future for a calibration of the mass-luminosity relationship. Even more dynamical masses will be required to extend the study 
of UCD physical properties to additional parameters, such as age, gravity, and metallicity. The study of UCDs would therefore 
greatly benefit from new high spatial resolution surveys dedicated to searching for new multiple systems, and from complementary 
monitoring programs targeting the shortest period binaries. 
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Center/California Institute of Technology, funded by the National Aeronautics and Space Administration and the National Science Foundation. 
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Table 2. Sample 



Name 


SpT 


I 




J 


H 


K 


Ref. 


2MASSW J0o50359+105715 


L6 






16.5 


1 J. 


.z 


1 zL 


.J 


Kirkpatrick et al. (1999) 


2MASSW J0920122-I-351742 


L6.5 


19, 


,4 


15.6 




7 


1 ^ 


Q 


Kirkpatrick et al. (2000J 


2MASSW J 1 146344+223052 


L2/L2" 






14.2 




■ Z 


1 

iZ, 


.0 


T7" 1 1- " 1 1-1 / 1 r^nrw 

Kirkpatrick et al. ( 1999) 


2MASSW J1426316+1557013 


M8/L1.5 


16. 


5 


12.9 


12. 


.2 


11, 


.7 


Oizis et al. (2000b) 


2MASSW J1311391+803222 


M7.5/M8 


16. 


,2 


12.8 


12. 


.1 


11, 


.7 


Gizis et al. (2000b) 


2MASSW J15344984-2952274 


T5.5 






14.9 


14. 


.9 


14, 


.8 


Bursasser et al. (2002) 


2MASSW J1728114+394859 


L7 






16.0 


14. 


.8 


13, 


.9 


Kirkpatrick et al. (20()0) 


2MASSW J233 1016-040619 


«L2 


16. 


,3 


12.9 


12. 


.3 


11, 


.9 


Gizis et al. (2000b_) 


2MASSW J2140293+162518 


M9 






12.9 


12. 


.3 


11, 


.8 


Gizis et al. (2000b3 


DENIS-P J035726.9-441730 


1^9/0.1.5" 


18. 


1 


14.6 


13. 


.5 


12, 


.9 


Martin etal. (1999J 


DENIS-P J100428.3-1 14648 


M9.5/L0.5'' 


18. 





14.9 


14. 


.1 


13, 


.7 


Martin etal. (1999) 


DENIS-P J144137.3-094559 


LI 


17. 


,3 


14.2 


13. 


.2 


12, 


.4 


Martmetal. fl999) 


2MASSW J 1047 127 +402644 


M8 






11.4 


10. 


.8 


10, 


.4 


Gizis et al. (2000a) 


2MASS1 J 1847034+552243 


M7 






1.9 


11. 


.2 


10 


.9 


Cruz et al. (2003) 



Note. — I,J,K magnitudes of the DENIS objects from the DENIS survey; I,H,K5 magnit udes of the 
2MASS and LHS objects from the 2MASS survey; I magnitudes for the 2MASS objects from lBouv et al.l 
l|2003); H magnitudes for the DENIS objects from the 2MASS survey. If not specified, the spectral type 
corresponds to that of the unresolved system. Unless specified, the spectral type(s) correspond to those 
given in the last column reference. 

"Spectral type from lMartm et al.l (120061) 



Table 3. Main characteristics of the instruments used in this study 



Instrument 


Filter 


Platescale 


Field of view 


A/D 


Resolution 






[mas/pixel] 


["] 


[mas] 


[mas] 


HST/WFPC2 PC 


F814W 


45.5 


44'.'2x44'.'2 


85 


60 


HST/ACS HRC 


F814W 


25" 


35'.'4x38'.'0 


85 


40 


HST/NICMOSl 


FUOM 


43.2 


15'.'7xl5'.'7 


115 


90 


HST/STIS 


LongPass 


50.8 


6'.'9x6'.'9 


75 


60 


VLT/NACO 


Ks 


13.3 


13'.'6xl3'.'6 


68 


40 


Gemini/Hokupa' a 


Ks 


20 


2a'5x20'.'5 


68 


50 


Subaru/CIAO 


Ks 


21.3 


21'.' 8x21'.' 8 


68 


60 


Lick/AO 


Ks 


76 


19'.'4xl9'.'4 


180 


110 



"For pipeline processed data with MultiDrizzle. The "raw" platescale of the ACS/HRC 
is 28x24.8 mas. 
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Dateof Obs. 


Instrument 


Sep. [mas] 


P.A ["] 


Amag 


Filter 


Ref." 


2MASSW J0850359-H105715 


01-02-2000 


HST/WFPC2 


157.2±2.8 


114.7+0.3 


1.47+0.09 


F814W 


(2) & (3) 


21-10-2002 


HST/ACS 


141.7+0.9 


124.6+0.36 


1.36+0.02 


F625W 


(1) 




HST/ACS 






1.21+0.02 


F775W 


(1) 




HST/ACS 






0.91+0.08 


F850LP 


(1) 


09-11-2003 


HST/NICMOS 


127.4+4.3 


129.0+1.8 


1.10+0.04 


FllOM 


(1) 


2MASSW J0920122+351742 


02-09-2000 


HST/WFPC2 


75.1+2.8 


248.5+1.2 


0.88+0.11 


F814W 


(2)&(3) 


19-10-2002 


HST/ACS 


<40 






F625W 


(1) 




HST/ACS 








F775W 


(1) 




HST/ACS 








F850LP 


(1) 


10-03-2003 


HST/STIS 


<60 






LongPass 


(8) 


22-03-2003 


VLT/NACO 


<60 






Ks 


(1) 


03-10-2005 


HST/ACS 


<40 






F814W 


(1) 


08-04^2006 


HST/ACS 


<40 






F814W 


(1) 


2MASSW'' J 1 146344+223052 


28-04-2000 


HST/WFPC2 


294.1+2.8 


199.5+0.3 


0.75+0.09 


F814W 


(2)&(3) 


08-06-2002 


HST/WFPC2 


284.8+2.8 


205.2+0.6 


0.53+0.09 


F814W 


(1) 


13-06-2002 


HST/WFPC2 


282.7+2.8 


205.0+0.6 


0.55+0.09 


F814W 


(1) 


05-05-2003 


HST/WFPC2 


280.5+2.8 


207.6+0.6 


0.55+0.09 


F814W 


(1) 


10-02-2003 


HST/STIS 


275.1+2.8 


205.5+0.6 




Longpass 


(8) 


13-11-2003 


HST/WFPC2 


276.5+2.8 


209.0+0.6 


0.56+0.09 


F814W 


(1) 


2MASSW J1426316+1557013 


20-06-2001 


Gemini/Hokupa' a 


152+6 


344.1+0.7 


0.78+0.05 


J 


(4) 




Gemini/Hokupa' a 






0.70+0.05 


H 


(4) 




Gemini/Hokupa' a 






0.65+0.10 


Ks 


(4) 




Gemini/Hokupa' a 






0.57+0.14 


K 


(4) 


19-07-2001 


HST/WFPC2 


155.6+1.7 


333.7+0.6 


1.40+0.09 


F814W 


(3) & (5) 




HST/WFPC2 






0.76+0.11 


F1042M 


(3) & (5) 


10-03-2003 


HST/ACS 


194.4+0.9 


341.9+0.3 


0.99+0.08 


F625W 


(1) 




HST/ACS 






1.22+0.08 


F775W 


(1) 




HST/ACS 






1.31+0.08 


F850LP 


(1) 


28-04-2003 


HST/STIS 


194.6+2.8 


341.6+0.8 




Longpass 


(8) 


22-06-2006 


VLT/NACO 


265.8+1.8 


342.9+0.8 


0.57+0.02 


Ks 


(1) 


2MASSWJ131 1391+803222 


30-07-2000 


HST/WFPC2 


300.4+3.9 


167.2+0.7 


0.39+0.07 


F814W 


(3) & (5) 




HST/WFPC2 






0.45+0.09 


F1042M 


(3) & (5) 


25-04-2002 


Gemini/Hokupa' a 


267+6 


168.15+0.48 


0.14+0.05 


K' 


(6) 


27-02-2003 


HST/STIS 


262.7+2.8 


170.4+0.6 




Longpass 


(8) 


2MASSW J15344984-2952274 


18-08-2000 


HST/WFPC2 


65+7 


1+9 


0.5+0.3 


F814W 


(7) 




HST/WFPC2 








F1042M 


(7) 


19-01-2006 


HST/ACS 


198.8+0.9 


15.0+0.1 


0.26+0.03 


F814W 


(1) 


11-04-2006 


HST/ACS 


190.7+0.9 


15.1+0.1 


0.31+0.03 


F814W 


(1) 


2MASSW J17281 14+394859 


12-08-2000 


HST/WFPC2 


131.3+2.8 


27.6+1.2 


0.66+0.11 


F814W 


(3)&(5) 


07-09-2003 


HST/NICMOS 


159.6+4.3 


66.8+1.8 


0.15+0.04 


FllOM 


(1) 


14-08-2005 


HST/ACS 


182.4+0.9 


82.9+0.3 


0.45+0.04 


F814W 


(1) 


18-05-2006 


HST/ACS 


188.7+0.9 


86.2+0.1 


0.59+0.03 


F814W 


(1) 


01-01-2006 


HST/ACS 


195.0+0.9 


88.6+0.1 


0.50+0.03 


1-814\V 


(1) 


DENIS-P J035726.9-441730 


21-04-2001 


HST/WFPC2 


97.5+3.9 


174.3+2.3 


1.23+0.11 


F675W 


(3) 




HST/WFPC2 






1.50+0.11 


F814W 


(3) 


21-08-2002 


HST/ACS 


103.9+0.9 


175.6+0.5 


1.09+0.02 


F625W 


(1) 




HST/ACS 






1.13+0.02 


F775W 


(1) 




HST/ACS 






1.14+0.02 


F850LP 


(1) 


03-01-2003 


HST/STIS 


103.9+2.8 


176.7+1.5 




Longpass 


(8) 


13-09-2005 


HST/ACS 


104.1+0.9 


175.5+0.5 


1.19+0.07 


F814W 


(1) 


31-05-2006 


HST/ACS 


91.5+5.4 


178.2+0.4 


1.11+0.04 


F814W 


(1) 



DENIS-P J100428.3-1 14648 
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Bouy et al.: Follow-up observations of binary ultra-cool dwarfs 
Table 4 — Continued 



Date of Obs. 


Instrument 


Sep. [mas] 


PA n 


Amag 


Filter 


Ref 




27-10-2000 


HSTAVFPC2 


146.0±3.9 


305.3+1.5 


0.25+0.07 


F675W 


(3) 




HSTAVFPC2 






0.66+0.11 


F814W 


(3) 


14-02-2003 


HST/STIS 


133.9±2.8 


315.2+1.2 




Longpass 


(8) 


DENIS-P J144137.3-094559 


16-01-2001 


HSTAVFPC2 


375.3±2.8 


290.4+0.4 


0.30+0.07 


F814W 


(3) 


22-05-2001 


HSTAVFPC2 


372.5±2.8 


291.3+0.4 


0.28+0.07 


F814W 


(3) 


20-01-2002 


HSTAVFPC2 


367.8±2.8 


292.5+0.4 


0.26+0.07 


F814W 


(1) 


29-03-2002 


HST/STIS 


367.8±2.8 


293.0+0.4 




Longpass 


(8) 


01-05-2002 


HSTAVFPC2 


365.0±2.8 


293.2+0.4 


0.27+0.07 


F814W 


(1) 


01-01-2003 


HSTAVFPC2 


362.6±2.8 


294.9+0.4 


0.27+0.07 


F814W 


(1) 


03-01-2004 


HSTAVFPC2 


355.6±2.8 


297.3+0.4 


0.26+0.07 


F814W 


(1) 


2MASSW J 1847034+552243 


10-07-2003 


Subaru/CIAO 


82±5 


91.1 + 1.4 


0.16+0.10 


Ks 


(9) 


03-04-2007 


Lick/AO+LGS 


170±7 


112.2+0.3 


0.27+0.15 


Ks 


(1) 


2MASSW J1047127+402644 


25-04-2002 


Gemini/Hokupa ' a 


122±8 


328.36+3.75 


0.50+0.15 


Ks 


(6) 










0.91+0.20 


H 


(6) 


03-04-2007 


Lick/AO+LGS 


106±14 


319.3+1.0 


0.6+0.4 


Ks 


(1) 










1.2+0.4 


H 


(1) 


2MASSW J2140293+162518 


20-09-2001 


Gemini/Hokupa ' a 


155±5 


134.30+0.5 


0.75+0.04 


K' 


(6) 


21-05-2001 


HST/WFPC2 


159.0±2.8 


132.4+0.3 


1.51+0.11 


F814W 


(3) & (5) 




HST/WFPC2 






1.38+0.11 


F1042M 


(3) & (5) 


27-06-2006 


VLT/NACO 


108.7±1.3 


205.7+1.6 


0.73+0.02 


Ks 


(1) 


2MASSW J2331016-040619 


06-05-2001 


HST/WFPC2 


577±2.8 


293.7+0.4 


3.90+0.17 


F814W 


(3) 




HST/WFPC2 






3.54+0.17 


F1042M 


(3) 


22-09-2001 


Gemini/Hokupa ' a 


573±8 


302.6+0.4 


2.44+0.03 


K' 


(6) 


20-06-2003 


VLT/NACO 


586.0+30 


290+3 




Ks 


(1) 



"Reference for the measureme nt: (1 ) T his work; (2) iReid et all (|2000|): (3) i Bouv etal] (|2003|): (4) 
Iciose et al. (2002); (5) Gizis et al. ( l2003h ; (ei lClose et all ( l2003h ; f7') lBurgasser et alj ( l2003h ; (8l lMartfn et all 
I2OO6) ; (9) ,Siegleretal.,(i2005il 



Note. — When several filters are available at the same epoch, the given separations and positions angle 
correspond to the average of the values measured in the different filters, and the uncertainties to the propa- 
gated uncertainties. 
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Fig. 4 Relative orbital motion of the multiple systems presented in this paper. The dates of the first and last epochs are indicated. 
The primary is represented with a large circle. Esimated uncertainties on the measurements are smaller than the symbols, unless 
specificied. The instrumental uncertainties are somclimes clearly dominating, as in the case of 2MASSW J1426316-I-1557013 
(the 2"'^ epoch deviates significantly, probably because of a large uncertainty in the P. A of the camera on the sky). The open 
triangle in 21VIASSW J233 1016-040619 and 21VIASSW J17281 14-1-394859 panels represents the position that an hypotetic unrelated 
background star would have had at the last epoch. 
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Table 5. Comparison of observed and proper motions 



Object Proper motion [mas/yr] 


Observed motion B/A [mss/yr] Ref. 


iiiVlrtoo VV JUO JUJ J7T iUJ / i J 


1 44 7+7 n 


12+4 (1) 


1104.71774-407644. 


701+4 


5+4 (2) 


7MA'\'\W Tl 1 4fi^444-77^nS7 


Qf. n+n s 


14+4 (1) 


2MASSW J1426316-I-1557013 


97±2 


23+2 (5) 


2MASSW Jl 3 1 1 39 1 +803222 


291±5 


16+5 (5) 


DENIS-P J144137.3-094559 


204+18 


16+4 (4) 


2MASSW J15344984-2952274 


268.8+1.9 


23+3 (2) 


2MASSW J17281 14+394859 


45.0+6.4 


33+2" (3) 


2MASS1 J 1847034+552243 


148+6 


26+3 (5) 


2MASSW J2331016-040619 


249+1 


235+33" (5) 


2MASSW J0920122+35 1742 




<6 


DENIS-P J035726.9-441730 




3+2 


DENIS-P J100428.3-1 14648 




12+6 


2MASSW J2140293+162518 




33+2 


"The amplitude of the proper motion and the observed motion are comparable, but the orientations are 


inconsistent. See also Fig.|4] 






Note. — Proper motions from (DlDahn et all (l2002h: f2')lTinnev et al.l (l2003h: f3)IVrba et al.l d2004l): (4) 


iJameson et alj 120071); (5) USNO-B.l catalog; Observed motions evaluated using Tablel4l using the most 


distant measurements together with the corresponding epochs, and assuming a linear motion. 


Table 6. Period estimates (in years) 


using different methods 


Object 


Kepler's Law 


Change in PA Change in sep. 




(at max. elong.) 


(circular face-on) (circular edge-on) 


2MASSW J0850359+105715 


38 


95 80 


2MASSW J1047127+402644 


11 


195 151 


2MASSW J 1 146344+223052 


68 


134 237 


2MASSW J1426316+1557013 


44 


1501 27 


2M AS SW J 1 3 1 1 39 1 + 803222 


61 


290 82 


DENIS-P J144137.3-094559 


120 


155 226 


2MASSW J 15344984-2952274 


16 


144 12 


2MASSWJ17281 14+394859 


21 


32 44 


2MASS1 J1847034+552243 


18 


64 14 


2MASSW J2331016-040619 


147 


15 220 


2MASSW J0920122+351742 


5 




DENIS-P J035726.9-441730 


8 


1319 260 


DENIS-P J100428.3-1 14648 


44 


84 111 


2MASSW J2140293+162518 


19 


24 64 
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